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Abstract 

We concern the problem of modifying the edge lengths of a tree in minimum total cost so that the 
prespecihed p vertices become the p-maxian with respect to the new edge lengths. This problem is called 
the inverse p-maxian problem on trees. Gassner proposed in 2008 an efficient combinatorial alogrithm to 
solve the inverse 1-maxian problem on trees. For the case p > 2, we claim that the problem can be reduced 
to finitely many inverse 2-maxian problems. We then develop algorithms to solve the inverse 2-maxian 
problem for various objective functions. The problem under b-norm can be formulated as a linear program 
and thus can be solved in polynomial time. Particularly, if the underlying tree is a star, the problem can 
be solved in linear time. We also develop 0(n log n) algorithms to solve the problems under Chebyshev 
norm and bottleneck Hamming distance, where n is the number of vertices of the tree. Finally, the problem 
under weighted sum Hamming distance is AP-hard. 

keywords: Location problem; p-Maxian; Inverse optimization; Tree. 


1 Introduction 

Location theory plays an important role in Operations Research due to its numerous applications. Here, we 
want to find optimal locations of new facilities. Location theory was intensively investigated, see Kariv and 
Hakimi [M [20], Hamacher Eiselt [12]. Recently, a new approach of location theory, the so-called inverse 
location problem, has been focused and become an interesting research topic. In the inverse setting, we aim 
to modify the parameters in minimum cost so that the prespecihed locations become optimal in the perturbed 
problem. Inverse median and inverse center problems are the popular topics in the held of inverse location 
theory. 

For inverse 1-median problems, Burkard et al. [^ solved the inverse 1-median problem on trees and the 
inverse 1-median problem on the plane with Manhattan norm in 0(n log n) time. Then Galavii [15] proposed a 
linear time algorithm for the inverse 1-median problem on trees. Burkard et al. |5] investigated and solved the 
inverse Fermat-Weber problem in O(ulogn) time if the given points are not colinear. Otherwise, the problem 
can be formulated as a convex program. For the inverse 1-median problem on a cycle, Burkard et al. |8] 
developed an O(n^) algorithm based on the concavity of the corresponding linear programmming constraints. 
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Additionally, the inverse p-median problem on networks with variable edge lengths is A^P-hard, see Bonab et 
al. [5]. However, the inverse 2-median problem on a tree can be solved in polynomial time. More particularly, if 
the underlying tree is a star, the corresponding problem is solvable in linear time. Sepasian and Rahbarnia |24j 
investigated the inverse 1-median problem on trees with both vertex weights and edge lengths variations. They 
proposed an 0(n log n) algorithm to solve that problem. While most recently papers concerned the inverse 
1-median problem under linear cost functions, Guan and Zhang [16] solved the inverse 1-median problem on 
trees under Chebyshev norm ang Hamming distance by binary search algorithm in linear time. 

Cai et al. m were the first who showed that although the 1-center problem on directed networks can 
be solved in polynomial time, the inverse problem is TVP-hard. Hence, it is interesting to study some special 
situations which are polynomially solvable. Alizadeh and Burkard [I] developed a combinatorial algorithm 
with complexity of O(n^) to solve the inverse 1-center problem on unweighted trees with variable edge lengths, 
provided that the edge lengths remain positive thoughout the modification. Dropping this condition, the 
problem can be solved in O(n^c) time where c is the compressed depth of the tree. For the corresponding 
uniform-cost problem, Alizadeh and Burkard [2] devised improved algorithms with running time 0(n log n) 
and O(cnlogn). Then Alizadeh et al. [3] use the AVL-tree structure to develop an 0(n log n) algorithm for 
solving the inverse 1-center problem on trees with edge length augmentation. Especially, the uniform-cost 
problem can be solved in linear time. For the inverse 1-center problem on a simple generalization of tree 
graphs, the so-called cactus graphs, Nguyen ang Chassein m showed the A^P-hardness. If the modification 
of vertex weights is taken into account, the inverse 1-center problem on trees can be solved in O(n^) time, see 
Nguyen and Anh [22] . Furthermore, Nguyen and Sepasian [28] solved the inverse 1-center problem on trees 
under Chebyshev norm and Hamming distance in O(nlogn) time in the case there exists no topology change. 
In principle, the problem is solvable in quadratic time. 

Although the inverse location problem was intensively studied, there is a limited number of papers related 
to the inverse obnoxious location problem. Alizadeh and Burkard |4] developed linear time algorithm to solve 
the inverse obnoxious I-center problem. Moreover, Gassner [13] investigated the inverse 1-maxian problem on 
trees with variable edge lengths and reduced the problem to a minimum cost circulation problem which can 
be solved in 0(n log n) time. 

We focus in this paper the inverse p-maxian problem on trees with p > 2. This paper is organized as 
follows. We briefly introduce the optimality criterion of the p-maxian problem on trees as well as formulate 
the problem under arbitrary cost function in Section [2] Section [3] concentrates on the problem under li-norm. 
We show that the problem can be formulated as finitely many linear programs. If the underlying tree is a star, 
we can solve the problem in linear time. We focus on Section |4] and Section [5] the problems under Ghebyshev 
norm and Hamming distance, respectively. We show that these problem can be solved in 0(n log n) time. 

2 The p-maxian problems on trees 

Given a graph G = {V,E), \V\ = n, each vertex vi G V associates with a non-negative weight Wi and each 
edge has a non-negative length Ig. The length of the shortest path P{u,v) connecting two vertices u and v 
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in G is the distance d{u, v) between these two vertices. A point on G is either a vertex or lies on an edge of 
the graph. As proposed by Burkard et al. [9], the p-maxian problem on G is to identify a set of p points, say 
X = {xi, X 2 , ■ ■ ■, Xp}, to maximize the maxian objective fnnction 


n 



By the vertex domination property, there exists a p-maxian of G that is the set of p vertices. Now we 
restrict ourselves in the case where the underlying graph G is a tree, say T. Burkard et al. [9] already showed 
that {a, b} is a 2-maxian of a tree T if P{a, h) is its longest path. We further show that this condition is also 
the necessary condition as follows. 

Theorem 1. (2-Maxian Criterion) 

Given two vertiees a,b on the tree T, then {a,b} is a 2-maxian of T if and only if P(a,b) is the longest path 
in the tree. 

Proof. For the sufficient condition, see [S]. Assume that P{a,b) is not the longest path of the tree, we will 
prove that {a,b} is not its 2-maxian as well. Let ruab be the midpoint of path P(a,b). By deleting an edge 
that contains rUab, we can deduce two parts L and R oi T which contain a and &, respectively. Then we get 

F{{a, b}) = J2veL Wvd{v, b) -f Wvd{v, a). 

Moreover, let P{s,t) be the longest path of the tree. We trivially get d{mab,a) < ma.x{d{mab,s),d{mab,t)} 
and diruabib) < max{(i(mab, s), (i(mob, t)}, and at least one of the two inequalities does not hold with equal¬ 
ity. Otherwise, it contradicts the assumption that P{s,t) is the longest path. Therefore we get d{v,a) < 
max{(i(z;, s),d{v, t)} for v G R and d{v, b) < max{(i(u, s), d{v, t)} for v G L, and at least one of the inequalities 
will not hold with equality. We can finally conclude that F{a, b) < F{s,t). □ 

We now reformulate the optimality criterion in Theorem [TJ Given two prespecified vertices a and &, and a 
leaf V, let Vab be the common vertex of three paths P{a,v), P{b,v), and P{a,b). The following result states 
the necessary and sufficient conditions for {a, b} to be the longest path in T. 

Corollary 1. (Longest path criterion) 

P(a, b) is the longest path ofT if and only if d{v, Vab) < d{a, Vab) and d(v, Vab) < d(b, Vab) for all leaves v in T. 

Next we formulate the inverse version of the p-maxian problem on trees. Given a tree T = (V, E) and a 
prespecified p-vertex. We can assume without loss of generality that the prespecified p vertices are the leaves 
of T. The length of each edge e in A can be increased or decreased by an amount pe or qe, i.e. the new length 
of e is = -^e +Pe — Qe and is assumed to be non-negative. We can state the inverse p-maxian on T as follows. 

1. The p-vertex becomes a p-maxian of the tree with respect to the new edge lengths i. 

2. The cost function C(p,q) is minimized. 

3. Modifications are feasible, i.e. 0 < Pe < Pe and 0 < qe < fe- 
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Burkard et al. [9] states that the set S oi p vertices is a p-maxian of T if and only if it contains a pair of 
vertices {a, b} such that P{a, b) is the longest path of T. Therefore, the inverse p-maxian problem on a tree 
can be reduced to many 2-maxian problems on this tree. From here on, we focus the 2-maxian problem on 
the tree T. 

Consider a monotone cost function C and {a, b} is a pair of leaves which are the prespecified vertices, we 
get the following property. 

Proposition 1. In the optimal solution of the inverse 2-maxian problem on T, it suffices to increase the 
lengths of edges in P{a, b) and reduce the lengths of edges in T\P{a^ b). 

By Proposition [U we can set Xe '■= Pe and Xe '■= pe for e € P{a,b), and Xe ■= Qe and Xe = pe for e G 
T\P{a, b). An edge e is said to be modified by an amount Xe if its modified length is set to £e '■= £e + sign{e)xe, 
where signie) = 1 if e € P{a, b) and sign(e) = —1 if e ^ P^a, b). Furthermore, denote by £ the set of leaves 
in the tree T and d by the distance measure with respect to new edge lengths i. In order to make P(a, b) 
the longest path of T, it must hold that d(v,Vab) < d{a,Vab) or d{v,Vab) < d(b,Vab) for all leaves v in T, see 
Corollary [T] After some elementary computations, we can formulate the inverse 1-maxian problem on T as 
follows. 

min C(a:) 

s.t. Xe+ Y >S{a,v), 

e^P{a,Vab) e^P{v,Vab) 

Y ^e+ Y ^S{b,v), 

eGP{b,Vab) e£Pi'u,Vab) 

0 < Xe < Xe, 

Here, 0(*, v) := d(v, Vab) — d(*, Vab), -*: = a,b, ior v G C are the gaps between distances. 

3 The problem under /i-norm 

Assume that modifying an edge e by a unit amount costs Ce, then the objective function in ([1]) under Zi-norm 
can be written as 


VvGC 

( 1 ) 

MvGC 

Ve e E. 


C(cc) — 'PjeCiE ^eXe- 

The inverse 2-maxian problem on trees under li-norm can be formulated as a linear program. It is therefore 
solvable in polynomial time. However, an efficient combinatorial algorithm is still unknown. 

We now consider the case where the underlying tree is a star graph with center vertex vq- We can directly 
deliver the following property from Corollary [TJ 

Corollary 2. (Optimality criterion) 

Given two leaf nodes a and b in a star graph S. Then, {a, b} is a 2-maxian of the star graph if and only if 
^ivo,a) cind £{vo,b) two largest edges. 
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For simplicity, we denote i{vo,a) and i(vo,b) by ia and and the cost to modify one unit length of (vo,a) 
and {vo,b) by Ca and Cb- By Corollary [21 we consider how to modify the length of egdes with minimum cost 
such that ia and it become the two largest edges in the star. Asume that ia < it, we analyze these two 
following cases. 

1. If fa S [ia,ib], we do not increase the length of it but decrease the length of e ^ (vg, b) with ie > ia to 
ia- We first presolve the problem as follows. 

• If ^ = min {ie — Xe} > ia, then increase ia by an amount ^ — ia- 

e^{vQ,h)-le>^a 

• If 6 = ia + X(^aa^a) < ie for e ^ {vo,b), then we decrease the length of ie by an amount ie — d- 

If the presolution is not possible, then the problem is infeasible. Otherwise, we can formulate the problem 
as a univariate optimization program. 

vain f {z) = Ca{z - ia) + ^ Ce{ie-z) 

e^{vo,b)-.ie>z 

where z = 4 and z e [4,4 + S(«o,a)]- 

2. If ia > ib, then the length of e where ie > ia must be decrease to 4 and the length of ib must be increase 
to ia- First, we have to increase ia an amount ib — ia and presove the problem as in the first case. If 
the presolution is not possible then the problem is infeasible. Otherwise, we can fomulate the problem 
as follows. 


min/(z) = (Ca + Cb)(z - 4) + ^ Ce(4--z) 

e:ie>z 

where z = ia = ib and z € [4,min{4 + 4 + S(«o,ft)}]- 

In both cases, the function /(z) is a convex function. According to Alizadeh and Burkard [3], the univariate 
optimization problem can be solve in linear time. To get the optimal solution of the inverse 2-maxian problem 
on S, we have to solve two problems and get the best one. Therefore, we attain the following result. 

Theorem 2. The inverse 2-maxian problem on a star graph ean he solved in linear time- 


4 Problem under Chebyshev norm 

We recall that modifying an edge e by one unit length costs Ce, the objective function in o under Chebyshev 
norm can be written as 


C{x) = maXe(zE{CeXe} - 

The solution method of the problem is based on the greedy modification. We define a valid modification 
of egde lengths in T with cost C as follows. 
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Definition 1. (Valid modification, see [23\ } 

A valid modification with cost C is to modify each edge of T as much as possible such that the cost is limited 
within C, i.e. we set 

,if CeXe > C 

,if CeXe < C. 

We can solve the problem by the following two phases. 

Phase 1: Find the interval that contains the optimal cost. 

We first sort the costs {ceXe\e^E in nondecreasing order and unite the similar values. Then we obtain a 
sequence of costs 



CiXi < C2X2 < ... < 

Here, m = 0{n). 

Then we aim to find the smallest index ig such that P{a, b) become the longest path by applying the valid 
modification with cost Ci^ Xi^. The index ig can be found by applying a binary search algorithm. If we apply 
the valid modification with cost CkXk but the path F(a, b) does not become the longest path then ig > k. 
Otherwise, we know ig < k. 

Let us analyze the complexity to find ig. In each iteration modifying the edge lengths of the tree costs 
linear time. Then we find the longest path of the tree in linear time (see Handler [15]), and compare the length 
of P{a, b) with the longest one in order to decide if Corollary [1] holds or not in linear time. As the binary 
search stops after O(logn) iterations, this procedure runs in 0(n log n) time. 

Assume that we have found the interval [cig-iXig-i,Ci,,Xi^] that contains the optimal cost. We apply a 
valid modification with cost dg-iXig-i and get the modified tree T. Also, we update the upper bounds of 
modifications. The next step is to define a parameter t G [0, and find the smallest value t such that 

P(a, b) becomes the longest path of T for a valid modification with cost afi. 

Phase 2: Find the minimizer t with respect to the optimal objective value in [0, 

After deleting all edges of the path P{a, b) we get a forest Consider a vertex v in the path P(a, b), 

V ^ a, b. A branch Ty is, by definition, a connected component of that contains the vertex v. Denote 
by C(Ty) the set of leaves in Ty except a and b. Now we try to modify the edge lengths of T so that there is 
no path P(v,v'), for v G P(a,b), v a, b and v' G C{Ty), is longer than min{(i(u, a), (i(u, 6)}. To solve this 
problem, we apply a valid modification with cost Cifi for t G [0,a;ip] and consider the modifying lengths of 
these following paths. 

• For a path P(v,v') with v' G CiTy) and d{v,v') > min{(i(w, a), (i(w, &)}, its modifying length is 
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d{v,v')-Ci„ . 

eGP{v,v'),Xe>0 ^ 

• For the paths P{v, a) and P{v, b), we obtain the modifying lengths 

d{v,a) + Ci^ ^ , 

eGP{v,a),Xe>0 ^ 

and 

d{v,b)-\-Ci^ ^ . 

eGP{v,b),Xe>0 ^ 


We first identify the value ti and ^2 such that 

ti := argminmax{(i(w,w') - Cio ^ , d{v, a) + ^ —}, (2) 

eGP{v,v'),Xe>0 ^ e(^P{v,a),Xe>0 ^ 

^2 := argminmax{d(z;,'!;') - Cio ^ —,d{v,b) + Ctg ^ —}, (3) 

eGP{v,v'),Xfi>0 ^ e^P{v,b),Xe>0 ^ 

where z; € P{a,b), v' € P(Tv), d(v,v') > min{d(v, a), d(v, b)} and t G [Q,Xi^]. As we have to find the minimizer 
of the upper envelop of linear functions, the algorithm in Sec. 5 of Gassner m will find ti and t 2 in linear 
time provided that all linear functions have been found. 

The smallest value t*, such that the modifying lengths of all paths P{v,v') for v G P{a,b)\{a,b} and 
v' G C{Ty) are not larger than min{(i(z;, a), (i(u, 6)}, is easily identified as t* := max{<i,t 2 }. 

We now analyze the complexity to find all linear functions in the expressions (I2|), ([3]). We number the 
vertices in the path P{a, b) from the left to the right side by a, vi, V 2 ,... ,Vk,b, where k = 0(n). Denote by 
E{S) the set of edges in a subtree S. We start from vertex vi, it costs 0{\E{Ty^)\) time (by a breath-first-search 
procedure) to identify all functions 

d(vi,v')-Ci„ V — 

^— Ce 

e^P{vi,v'),Xe>0 


for v' G C{Tvi) and 0{\E{P{a,b))\) time to identify 

d{Vl,-k) + Cig 


e&P(vi ,*),Xe>0 

for -k = a,b. 

Next we have to find all linear functions with respect to vertex V 2 ■ To identify the functions 

d{v2,-k)+Ci„ ^ 

eGP{v2.,*),Xe>0 


t 

Cf> 


for = a, 6, we have to exchange the modification edge (ni, ^ 2 ) that contribute to the augmentation of path 
P{v 2 ,a). We continue the process for vertex V 3 ,V 4 ,... ,Vk similarly. 

In conclusion, it costs 0{J2vGP{a b)\{a b} + 2|£’(P(a, b))\) = 0{n) time to complete the procedure 

to find all linear functions in ([2]) and (jS]). 


Theorem 3. The inverse 2-maxian problem on a tree under Chebyshev norm can be solved in 0{nlogn) time. 
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5 Problem under Hamming distance 


First, let us focus on the problem under bottleneck Hamming distance. Assume that modifying one unit length 
of edge e costs Ce, we aim to minimize the following objective function 

maxegs{ceH(a;e)}. 

Here, iJ is a Hamming distance and defined by 

f 0 , if 0 = 0 

H{9) := 

I 1 , otherwise 

By the special structure of the Hamming distance, the objective function receives finitely many values, say 
{ce : e S E}. Therefore, we can solve the problem by finding the smallest value in {ce '■ e € E} such that the 
optimality criterion in Lemma [I] holds. 

Number the edges in T by 1,..., m, and the corresponding costs are ci,..., Cm, for m = n — 1. Let us first 
sort the costs {ce : e G E} increasingly and get without loss of generality a sequence 

Cl < C2 <■■■ < Cm • 

Now we apply a binary search algorithm to find the optimal cost. We start with the cost Cfe, k = 

We modify all the edges l,2,...,fc — 1. If the optimality criterion holds, we know that the optimal value is 
less than or equal to Ck- Otherwise, it is larger than Cfe. In each iteration we recomputing the length of the 
longest path in linear time (see Handler [18]) and compare it with the length of P{a,b). Moreover, as the 
binary search stop after O(logm) iterations. Phase 1 runs in 0(m log m) = 0(n log n) time. 

Theorem 4. The inverse 2-maxian problem on trees under bottleneck Hamming distance can be solved in 
0{n log n) time. 

For the problem under weighted sum Hamming distance, the objective function can be written as (xe). 

We can easily reduced the Knapsack problem into an inverse p-maxian problem under weighted sum Hamming 
distance in polynomial. Therefore, we get the following result. 

Theorem 5. The inverse p-maxian problem on tree under weighted sum Hamming distance is NP-hard. 

6 Conslusion 

We have addressed the inverse p-maxian problem, p > 2, under various objective functions. It is shown that 
the problem can be reduced to p^ many 2-maxian problems. Then we have formulated the inverse 2-maxian 
problem on trees under li-norm as a linear program and solved the problem on star graphs in linear time. 
Furthermore, the inverse 2-maxian problem under Chebyshev norm and Hamming distance is solvable in 
0(n log n) time, where n is the number of vertices in the tree. For furture research topics, we will consider the 
inverse maxian problem on other classes of graphs, e.g., cacti, interval graphs, block graphs, etc. 
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